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' Abstract. We make a direct connection between the construction of three dimen- 

sional topological state sums from tensor categories and three dimensional quantum 

' g rav hy by noting that the discrete version of the Wheeler-DeWitt equation is ex- 

actly the pentagon for the associator of the tensor category, the Biedenharn-Elliott 
identity. A crucial role is played by an asymptotic formula relating 6j-symbols to 
rotation matrices given by Edmonds. 

m : 
o ■ 
OS ■ 

O | Introduction 
SO . 

It has been known since the work of Ponzano and Regge [1] that the quantum 
theory of gravity in three Euclidean dimensions can be described by means of the 
evaluation of spin networks, which are combinations of representations of the Lie 
group SU(2) under the tensor product in the category of representations of SU(2). 
J> ! The ideas were elaborated by Hasslacher and Perry [2]. 

The evaluation of a spin network on a sphere can be calculated by triangulating 
the three-dimensional ball bounded by the diagram and labelling all interior edges 
with representations. Ponzano and Regge gave a formula as a sum over all such 
labellings of the product of the 6 J symbols associated to the tetrahedra of the 
triangulation. The edge labels j are identified with half-integers. 



\& (boundary edge labels) = 

E II 6j-symbol) J] (2j + 1) (1) 

interior edge labels Vtetrahedra / interior edges 

where the integer x is a function of the edge labels. If the sum is finite, then it 
is independent of the triangulation in the interior of the ball. A more detailed 
account of this has appeared in [3]. The formula can also be applied to obtain a 
wavefunction for other 3-manifolds more general than the ball. 
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Subsequent developments have put this formula into a broader context. Turaev 
and Viro [4] showed that the analogous formula, using the truncated category of 
representations of a quantum group at a root of unity gives rise to a three di- 
mensional topological quantum field theory (TQFT), closely related to the Chern 
Simons theory of Witten [5] . 

Barrett and Westbury [6,7], have demonstrated that the natural generalisation 
of this construction is a tensor category with unit and duals and only a very weak 
and natural assumption on the left and right traces. In particular the braiding or 
symmetry of the tensor product is unnecessary. The useful generalization of 6J 
symbols is the associator isomorphism of the tensor category. 

Thus, we have what seems to be a rather mysterious connection between a very 
abstract branch of mathematics, namely the theory of tensor categories, and the 
very deep physical problem of quantising general relativity. 

The purpose of this paper is to make this connection more direct, and possibly 
somewhat clearer. What we show is that the Wheeler-DeWitt equation for 3d 
general relativity reduces directly to the pentagon relation for the tensor category 
we are using, in the presence of an asymptotic formula for thin tetrahedra. Since 
the pentagon [8] is just the coherence condition for the associator of the tensor 
product of the category, we are identifying the laws of motion of general relativity 
directly with the most fundamental part of the theory of tensor categories. 

This allows us to focus our efforts to make physical interpretations of theories 
based on other tensor categories on the discovery of the thin tetrahedra asymptotic 
forms of the associator isomorphisms. 

This paper should be thought of as shedding some light on the algebraic pro- 
gram for quantising general relativity outlined in [9]. The great desideratum, of 
course, would be to find an analog of the categorical formula which worked in four 
dimensions. We have not yet worked out the four dimensional argument analogous 
to the one given in this paper but we believe it should be possible to do so, and that 
it might help to recognize whether an algebraic construction of a four dimensional 
state sum is related to general relativity. 

The Ponzano-Regge wavefunction is defined for triangulations of a 3-manifold 
with boundary where the state-sum (1) is a finite sum. The wavefunction is then the 
Turaev- Viro partition function with the parameter q set to 1, up to a normalisation. 
This shows that when the Ponzano-Regge wavefunction is defined it depends only 
on the boundary data and the topology of the 3-manifold, being independent of 
the particular triangulation. However the Ponzano-Regge formula does not define 
a TQFT because for an arbitrary triangulation the formula (1) may be infinite. 

In any TQFT defined by a state sum in a similar manner to the Turaev- Viro 
state sum, there is a vector space V determined by a triangulation of a surface. 
The partition functions, or wavefunctions, all lie in a particular subspace S C V. 
While the space V depends on the triangulation, the subspace S is independent of 
the triangulation. 

One can see on general grounds that S is determined by a set of projectors n v , 
one associated to each vertex v in the surface, and that S is determined by the 
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equations 

TYylf) = if) (2) 

being simultaneously satisfied. The equations are obviously candidates for the 
analogues in this theory of the Hamiltonian constraint equations which are often 
supposed to characterise diffeomorphism invariant quantum field theories. 

In this paper, we present constraint equations in the form of (2) for the Ponzano- 
Regge wavefunction. However the argument is independent of the general formalism 
of TQFT. We give an argument that the equations are a quantisation of the ap- 
propriate discrete version of the Wheeler-DeWitt equations, thus identifying the 
Ponzano-Regge wavefunction as belonging to a quantisation of a specific action, 
namely Einstein gravity. 

THE CONSTRAINT EQUATIONS 

The Einstein equation for a metric in three dimensions is that it is locally flat. 
On the boundary surface E of the manifold, this implies that there are constraints 
between the extrinsic curvature K ab and the metric g ab of E. Assuming that the 
three-dimensional metric is positive definite, these are 

V a K a b - V b K = 

and 

K 2 - K ab K ab -R = 0, (3) 

where K = K ab g ab , V is the connection for g, and R its scalar curvature. The 
conventions are those of [10]. 

The Einstein action can be used to cast these equations in the framework of 
Hamiltonian mechanics [11], using the momentum density 

iv ab = (Kg ab - K ab ) Q (4) 

where O is the metric volume measure. Then the constraints are 

(tt 2 - n ab n ab ) IT 1 - RQ = (5) 

and 

V a yr a 6 = 0. (6) 

For each tangent h ab to the space of metrics on E, canonical quantisation associates 
to the momentum 

<t>= f n ab hab (7) 

the first order differential operator corresponding to h ab . 

Making this substitution in the Hamiltonians determined by (5) and (6) does not 
generally make sense in this infinite-dimensional configuration space. However this 
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system can be quantised by changing variables and fixing the gauge in a particular 
way [12]. 

Alternatively, one can make progress by replacing the space of metrics on E by a 
finite-dimensional approximating space. Then canonical quantisation gives partial 
differential equations, with some ambiguity due to the different possible orderings 
of the operators. 

Our quantisation procedure is to use the finite-dimensional space of metrics de- 
termined by Regge calculus [13] and regard quantisation as giving an approximation 
to the equations satisfied by the wavefunction, which is already defined. Thus we 
are not using quantisation to define a quantum theory, rather to compare a quantum 
theory with a classical mechanical one. 

The momentum cf> survives the discretisation procedure, and for a metric flucta- 
tion given by changing one particular edge length I, is just the parameter for the 
extrinsic curvature at the edge discussed by Hartle and Sorkin [14] . This identity is 
specific to three dimensions, the momentum and the extrinsic curvature being log- 
ically distinct quantities. If E is the boundary of a 3-manifold M with a Euclidean 
signature metric, this is the angle of rotation of the outward unit normal from the 
triangle at one side of an edge to the normal to the triangle at the other side. It is 
positive if the edge is convex and negative if the edge is concave. 



b 




Figure 1. Part of E 



Consider a vertex in E which belongs to exactly three triangles, such as appears 
on the boundary of a tetrahedron (figure 1). The connection in M gives a map 
TM P — > TM p f if a path from p to p' in E is given which avoids vertices. If an 
orthonormal frame is chosen at p and p' , then this map determines a rotation in 
SO(3) 



rotation 



R 3 



frame at p 



frame at p 

-> TM P , 
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Consider a path which circulates the vertex via six points at which canonical 
frames are chosen, as shown in figure 2. At each point, vector z is chosen to lie 
along the adjacent edge tangent to E and x is chosen rotated by a local orientation 
of the surface. The frame is completed by vector y which is the outward normal to 
E. The six rotation matrices are Ri 2 , -R23, -R31, which are rotations through #12, 
6*23, #31 about the y-axis, and Ki, K 2 , which are rotations by 0i, 02, 03 about 
the z-axis. The i?'s can be described as the rotations which take one edge into the 
next and the K } s are the rotations which make one face parallel to the next. 

As for the smooth case, the Einstein equations in three dimensions are that the 
metric is flat. The condition that M is flat is that the holonomy in M is the identity, 

h = R 2 ^K 3 R 31 K 1 R 12 K 2 =Ie SO(3). (8) 




Figure 2. Path in E and choice of frames 



This gives three equations per vertex of E which constrain the metric of E and 
its extrinsic curvature. These are the key equations which we wish to quantise, to 
give our discrete version of the Wheeler-DeWitt equation. 

For the rest of this section there is an argument to demonstrate that the equations 
(8) are the discrete analogues of the Hamiltonian constraint equations. However 
the development of the rest of the paper is independent of this argument. 

This can be seen by considering small perturbations of the geometry around a 
flat configuration 0i = 02 = 03 = 0. One gets three equations, corresponding to the 
three standard basis elements X, Y, Z of the Lie algebra of SO (3), which generate 
rotations about the x, y and z axes. At the flat configuration, the two-dimensional 
holonomy is the identity, 

R23R31R12 = I e SO(3), 

and a small perturbation results in a small deficit angle 5. Each edge has associated 
a Lie algebra element 

L x = (cos6> 12 )Z+ (sin6>i 2 )X 
L 2 = Z 

L 3 = (cos 6 23 )Z- (sme 23 )X 
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Expanding (8) for small 0i, 02, 03, and keeping only the lowest-order non-zero 
terms, one gets 

4>\Li + 2 2v2 + 03-^3 = 

and 

<5 = 7^0102 sin ^12- 

The first of these is two equations, linear in the momenta, given by the coefficients 
of X and Z. It is a divergence equation, and corresponds to (6). The second 
equation, quadratic in the momenta, corresponds to (5), and is symmetrical in the 
three edges by virtue of the first equation, giving the area of a triangle with edges 
01) 02, 03- It is the coefficient of Y on expanding (8). 

Quantisation 

The Ponzano-Regge wavefunction is a function of the edge labels on the boundary 
E of M. The length of an edge is defined to be the value of the label plus one 
half. Thus the edge lengths are discrete, taking the values j • + \ : where j G |Z, 
j > 0, indexes the representations of SU(2). These lengths determine a Euclidean 
signature metric for E. We consider the equations satisfied when the edge lengths 
Zi, Z2, I3 around a vertex at which the three edges meet are varied, the other edges 
in E having fixed lengths. The simplest example is for a single tetrahedron, when 

^{hMM) = { a h I (9) 

the 6j-symbol, for fixed a,b,c. 




Figure 3. Crossing move 



Actually the case where three edges meet at a vertex suffices to handle the general 
the cases where more than three edges meet can be reduced to this by the 
crossing move, shown in figure 3, which is an invertible linear map V — > V and so 
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generates no new equations [3]. Also, any vertex at which three edges meet reduces 
to the case (9) because the Ponzano-Regge wavefunction factorises as 

il>(h, h, h) = | ^ i 2 | ^(a, b, c, other edges) 

by the excision move, shown in figure 4 [3]. Alternatively, this can be seen as an 
expression of the invariance of the wavefunction under a change of triangulation of 
the interior of M. 



b b 




Figure 4. Excision move 



Our result for the Ponzano-Regge wavefunction is that for each representation 
J it satisfies an equation 

J 

D%%^(hMM) = (-l? J £ DiXD^^ih + muh + m^h + mz) (10) 

m\ = — J 

In this formula, there are three square matrices of numerical coefficients D mn , of 
dimension 2J + 1. These matrices depend on the edge lengths (Ii,l2,l3,a,b,c). 
More specifically, Z)( 23 ) depends on (I2, h, a), Z)( 31 ) depends on Is, 6), and Z^ 12 ) 
depends on (/i,/ 2 ,c), so that each matrix is associated to one of the triangles in 
Figure 1. 

Proof of formula (10). Apply the Biedenharn-Elliott relation to the complex shown 
in figure 5. This gives 

\ a T r i 3 , i 2 U 7 a 1 7 C |= E (-lTO + 2m 1 + l) 

x mi = -J 

[b h l 3 1 f c l 2 h 1 f a 6 c 1 

\ J Z 3 + m 3 Zi + mi J \ J li + mi Z 2 + "^2 J \ Zi + mi / 2 + m 2 / 3 + m 3 J 
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Figure 5. Complex 

with x = a + b + c + h + h + h + J + h + mi + h + ^2 + h + "7-3 • This is equation 
(10), with 



DgfLs = (-l)' 2+i3+a+J+m3 (2/2+2m 2 +l)^(2/ 3 +2m 3 +l)^ | 



1 I a I3 I2 
J l 2 + m 2 I3 + m 3 



(and likewise for L>( 31 ), D^). 

This equation can now be compared to (8). We claim that it is a quantisation 
of (8). First, (8) expressed in any irreducible representation J. It can be written 
very explicitly using the weight basis in which the ICs are diagonal. 



K n = 



\ 



iJ<t>r, 



,*(J-1) 



\ 



e —iJ(f> n J 



and the i?'s then have the real matrix elements denoted d J mm ,{8) in the angular 
momentum literature, related to the Jacobi polynomials. For example, 

A(a\- ( cos6 / 2 *™e/2\ 

[0) \-sm6/2 cos6»/2 J ' 

/±(l + cos6>) ^sin6> ±(l-cos0)\ 
— 4= sin 9 cos 9 sin 9 



\\{\-cos9) -^sinfl ±(l + cos0)/ 



d\9) = 
Equation (8) becomes 

Cms (*23)= £ e^d^^e^di^iO^e^ 



(11) 
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According to the formula given by Edmonds [15], -Dm 2 m 3 is asymptotic to 

d J m 2 m 3 (023) = ^ 3m2 (-023) 

for fixed m2,W3, J but a,l2,h — > °o. Similar formulae hold for D^ 12 \ D^ 31 K For 
this approximation to be valid, it is therefore necessary to assume J is much smaller 
than the six edge labels appearing in Figure 1. 

Then the equation (10) is a quantisation of (11) in which e im ^ is replaced by the 
shift operator T 

Til>(l) = i/>(l + m), 

which was suggested as a heuristic by Ponzano and Regge. 

The other noteworthy feature is the factor (— 1) 2J . This is —1 for odd spin 
representations and gives the extension of (8) to SU(2) as the equation 

h sv ^ = -I. 

This is the natural extension as the product of the matrices in (8) gives a rotation 
which is exactly one full turn, which lifts to —I in SU(2). This phenomenon is 
discussed in [16]. 

Remarks. Our quantisation of the discrete Wheeler-DeWitt equation gives an 
appealing physical explanation of the Ponzano-Regge wavefunction by an analogy 
with other quantised systems such as particle mechanics. But it is not clear exactly 
what the mathematical 'point' of our work is; can anything now be proved that was 
not known before? 

The connection between spin networks and 2+1 dimensional gravity has al- 
ready been made in [1,2]. However this proceded in an indirect manner from the 
Biedenharn-Elliott relation to form a 3-term recursion relation for 6j-symbols, and 
then to a semi-classical limit. We have shown here that in the Hamiltonian picture 
the analog of the Wheeler-DeWitt equation is exactly the Biedenharn-Elliott rela- 
tion itself. This has two advantages: firstly it makes a more direct connection with 
the ideas of topological quantum field theory, since it focuses on the consideration 
of a manifold with boundary, and secondly it suggests an easier and more direct 
approach to making physical interpretations of other topological state sums. We 
note that the recursion relation occurs as a special case of our formula, namely 
when ni2 = = and J = 1. 

Our quantisation appears to give more information than the semiclassical formula 
for the 6j symbols given by Ponzano and Regge. For example, it gives an expla- 
nation for the discreteness of the edge length labels, which is not required by the 
semiclassical formula. Normally one would quantise the variable conjugate to the 
length / by <fi — > d/dZ. Then / would have to have a continuous spectrum. However, 
we are interested in quantising exp(im</>), which, making the above substitution, 
gives the shift operator, given by the action I — > I + m. As m is a half-integer, this 
operator natually acts on functions defined only on the discrete set of points |Z. 

Although we have considered quantising the formulae obtained from 3-dimen- 
sional metrics with a Euclidean signature, one can also start with the corresponding 
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Lorentzian formulae. Then there are some sign differences in (3) to (5), which are 
reflected in the fact that (8) becomes an equation in SO(2, 1) with the -fT's boosts 
and the i?'s unchanged, c.f. [17]. The quantisation replaces the matrix element e J< ^ 
with a shift operator, and one obtains the same asymptotic approximation to (10). 
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